Superfluidity versus Anderson localization in a dilute Bose gas 
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We consider the motion of a quasi one dimensional beam of Bose-Einstein condensed particles in a 
disordered region of finite extent. Interaction effects lead to the appearance of two distinct regions of 
stationary flow. One is subsonic and corresponds to superfluid motion. The other one is supersonic, 
dissipative and shows Anderson localization. We compute analytically the interaction-dependent 
localization length. We also explain the disappearance of the supersonic stationary flow for large 
disordered samples. 
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Interference effects have a strong influence on the 
transport properties of phase coherent systems. In par- 
ticular in one dimension (ID) in presence of disorder, 
they lead to Anderson localization (AL), revealed for in- 
stance by a transmission decreasing exponentially with 
system size. On the other hand, interactions in a low 
temperature Bose system may lead to superfluidity (SF), 
i.e., perfect transmission. Understanding the interplay 
between these two contrasting tendencies is one of the 
key issues of the physics of interacting disordered Bose 
systems. Recent advances in creating and manipulat- 
ing guided cold atomic vapors and atom lasers (see, e.g., 
Refs. [H, 0]) offers new prospects for studying such trans- 
port phenomena in Bose-Einstein condensates (BEC). 
BEC systems are of particular interest in that respect 
because they are almost perfectly phase coherent, and 
interaction between their components can be easily mod- 
ified and modeled. Recently, great experimental efforts 
have been devoted to the identification of AL in these 
systems 0, 0, 0|- Though not yet observed, the possibil- 
ity of AL has been pointed out theoretically in the non 
interacting regime [y, 0] , while interactions may play a 
role in the localization scenario Q. 

In this Letter we study a quasi fD, weakly interact- 
ing BEC, propagating through a disordered potential. In 
this context, localization has been theoretically studied 
mainly for effective attractive interactions (see, e.g., Q 
and references therein), with less attention on the re- 
vulsive case that we consider here (see, however, Refs. 

[13, EH)- I n the regime of weak disorder and weak 
velocity, the condensate depletion and the destruction of 
SF have been considered theoretically in Refs. [T^, while 
the strong disorder regime has been analyzed numerically 
in Ref. 

Our main result is a new, global picture of coherent 
transport of an interacting BEC through disorder, with 
a clear characterization of the different transport regimes 
where SF and AL can be observed. Our findings are 
summarized in Fig. [TJ The different types of flow are 
displayed as a function of v = V/c and L/£, where V is 
the condensate velocity relative to a disordered potential 
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FIG. f: (color online) Transport of a quasi ID BEC with 
normalized velocity v through a disordered potential of length 
L [Eq. ([2|, A = rii£ = 0.5]. Dark region: time dependent flow; 
light blue/gray regions: stationary flow. Note the enlarged 
scale for v € [0,1]. 



of spatial extension L, c is the speed of sound, and £ the 
healing length. At low velocities (subsonic regime), the 
flow is SF, as expected on the basis of Landau's criterion. 
The density profile is stationary, but locally modified by 
the external potential at the expense of the creation of 
a normal fluid fraction, that we compute perturbatively, 
extending previous results obtained in the limit v — > 
[l2| to finite velocities [Eq. |5])]. In the opposite super- 
sonic regime, a region of stationary flows also exists, but 
in this case energy dissipation occurs. In this domain, 
depending on the extent of the disordered potential, the 
system is either in an ohmic or in an AL regime (respec- 
tively characterized by a transmission decreasing linearly 
or exponentially with L, see below). In between the high 
and low velocity stationary regimes there is an interme- 
diate domain, centered around v ~ 1, where the flow is 
time dependent. 

Consider a BEC at rest in a 3D guide where a disor- 
dered potential U of finite spatial extension L is moved in 
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the longitudinal x direction at constant velocity V > 0. 
The transverse confinement is such that the "ID mean 
field regime" [T3| holds. The system is then described 
by an order parameter of the form exp{— iut/h}ip(x, t), 
which is solution of 



where 



. dip h 2 d 2 ip 
dt 2m dx 2 



[U(x-Vt)+g\1>\*-n]l>. (1) 



In the case of particles experiencing an effective repul- 
sive interaction characterized by the 3D s-wave scattering 
length a > 0, the constant g reads g — 2hv±a, where v\_ is 
the frequency of the transverse harmonic trap [l5| . In the 
stationary regime, where the flow is time-independent in 
the frame moving with the potential, tp depends on x 
and t only through the variable X = x — Vt. The ap- 
propriate boundary condition is ip(X — » — oo) = y/no 
(where no is a constant) [16|. The condensate is char- 
acterized by a chemical potential u = gno, a speed of 
sound c = (fi/m) 1 / 2 and a healing length £ = h/(mc). 
We consider a regime where the typical value of the dis- 
ordered potential is much smaller than /i + mV 2 /2. This 
regime is easily reached experimentally and is important 
for our purpose because it corresponds to a range of pa- 
rameters where AL is not blurred by effects connected to 
"fragmentation of the condensate" (see, e.g., Ref. and 
references therein). 

The picture we draw in the present Letter is generic 
and does not depend on the specific structure of the dis- 
ordered potential U. We illustrate our findings with a 
potential of the form 
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U(x) = A fi g S(x - x n ) 



(2) 



n=l 



where the the random, uncorrelated positions of 

N{ delta scatterers uniformly distributed over a length 
L (with a density n; = Ni/L). A > is the dimension- 
less strength of a scatterer. This potential has a mean 
value (U) = Xn^m andacumulant (U(x)U(0)) — (U) 2 = 
(h 2 /m) 2 a 6(x), where a = A 2 n\ £ -2 . From what is known 
in the case of Schrodinger equation, this potential is typ- 
ical insofar as localization properties are concerned. 

Numerically, we solve Eq. (TTJ) using the potential ^ 
for given u, to, m, A and V. We consider 100 realizations 
of the random potential. The fraction / of potentials 
for which a stationary solution exists is plotted in Fig. Q] 
using a gray scale (dark, / = 0; light blue/gray, / = 1) as 
a function of the normalized variables L/£ and v = V/c. 
This normalization rescales all interaction effects. In the 
following we characterize the different stationary regimes 
represented in Fig. Q] 

Let us start with the subsonic stationary regime. 
There, the density profile is easily computed perturba- 
tively to be of the form n{X) = no + 5n(X), with [3, 17| 



Sn(X) ~ - 



^ / dyU(y)exp{-2K\X-y\} 



(3) 



K= -| c2 _y 2 | 1/2 
h 1 1 



(4) 



is an effective wave vector. From Eq. ([3]) Sn(X — * ±oo) = 
and the density is only perturbed in the region of the 
potential. This corresponds to perfect transmission of 
the condensate through the disordered potential. There 
is no energy dissipation nor drag exerted on the potential, 
a characteristic feature of SF [1S|] . We have verified these 
properties numerically and also checked the stability of 
the SF flow. There is also a non superfluid fraction, since 
a momentum P = hj dx5nd x S [19| can be associated 
to the flow (where S is the phase of the order parame- 
ter). In the spirit of Landau's approach for determining 
the normal fraction in liquid 4 He (see also (13)) one can 
associate to this momentum a normal mass M n — P/V. 
In the perturbative limit where Eq. |3j) holds one gets 
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2H 4 k 3 L J R 2 



dyidy 2 U(y 1 )U(y 2 )K(y 2 -y 1 ) , (5) 



where K{y) = (1 + 2k|j/|) exp{— 2k|j/|} and M — mn L 
is the mass of a region of size L of the unperturbed con- 
densate. The ratio M n /M corresponds to the normal 
fraction present in the disordered region. For a potential 
of type (|2|) , at zero velocity, in the scarce impurities limit 
(1 > ni£), Eq. © yields M n (0)/M = \ 2 n i $,/2 1 in accor- 
dance with the findings of Refs. [HI, The present 
approach allows to extend this result to finite velocities, 
yielding (for k > n { ) : M a (v)/M n (0) = (1 - v 2 )- 3 / 2 . 

The situation changes drastically as the velocity in- 
creases. Starting from the SF subsonic region (with fixed 
U and L), at some critical velocity the system enters a 
domain where the flow becomes unsteady (dark region 
in Fig. [T]) . The velocity at which this transition occurs 
depends very much on the specific form of the potential 
U (see, e.g., the discussion of the influence of the value 
of A for a regular array of equally spaced impurities in 
[2ll]). If v further increases one gets into a new station- 
ary region (now supersonic). Within this region, if v is 
large enough, one reaches a domain (qualified as ohmic 
in Fig. [T]) where perturbation theory holds, yielding [l6| 



6n(X) 



2 m no 
I 2 



^ f dyU(y) sm[2rt(X - y)] 

K J-oc 



(0) 



where k is given by Eq. (gj). In the region upstream 
with respect to the disordered potential (X — > +00), 
Eq. (|6|) corresponds to an oscillatory wave, i.e., to a non- 
perfect transmission, indicating the dissipative nature of 
the flow. Equation ([6]) allows to express the disorder av- 
eraged transmission coefficient as [lOj : 



(T)~l 



L 



where L\ oc (n) = — . (7) 
a 
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FIG. 2: Statistical properties of the transmission T at v = 30 
for potential ([2| with A = ni£ = 0.5 [except plot (a) where 
v = 3 and A = 0.1]. Dashed line: analytical result; solid 
line: numerics. Top row: averaged transmission. Analytical 
and numerical results are almost indistinguishable. The non 
interacting results are also represented in plot (a). Bottom 
row: Probability distribution of T. Case (c): ohmic regime, 
L/L\ oc ~ 0.1; case (d): AL regime, L/L\ oc ~ 2.4. 



The probability distribution of T can be shown to be 



P(T) = %£exp 



(1 _ T) iloc 
Lj 



(8) 



The validity of Eqs. (|7l8p is confirmed by comparison 
with numerical results in Figs.[Ha) and^c). These equa- 
tions hold for potentials of the form @, for Gaussian 
white noise potentials as well as for correlated noise, but 
in the latter case the expression of L\ oc (k) is modified. 
In this regime the transmission decreases linearly with 
L, this is why we call it ohmic. Eqs. (|6|7|8p are valid 
in the perturbative regime \Sn\ <C no, or equivalently 
L <C L\ oc (k). Starting from the ohmic regime, keeping 
fixed v > 1 and increasing the length L of the disordered 
potential, one enters the interesting non-perturbative do- 
main L > L\ oc where AL (i.e., an exponential suppression 
of the transmission with disorder length) is expected. In 
this regime, we will show [cf. Eq. (fTT|) below] that the 
quantity L\ oc appearing in Eq. ([7]) corresponds indeed to 
the localization length. 

We now derive a non-perturbative treatment of the 
supersonic stationary flow, valid in the AL regime. We 
consider a disordered potential of type @ and look for 
stationary solutions. Between two impurities (x n and 
x n+ i say) the random potential is zero and |TJ) admits a 
first integral of the form [l(| 



i 2 

2 



V(^) - n-[.K>Y)] = r; 



1 ! 



(9) 



where £y is a constant, A — \i/j\/^/no and W(A) = 
\{A 2 - 1)(1 + v 2 - A 2 - v 2 /A 2 ). Eq. © has a Hamilto- 
nian form, expressing energy conservation for a fictitious 



FIG. 3: Upper panel: W as a function of A = \tj>\/y/rio (drawn 
for v = 4). A (= 1) and Ai are the zeros of dW/dA The 
fictitious particle is initially at rest with E°i = 0. The value 
of E c \ changes at each impurity. The lower panel displays 
the corresponding oscillations of A(X), with two impurities 
(vertical dashed lines) at x\ — and X2 = 4.7 £. 



classical particle with "mass" £ 2 , "position" A, "time" 
X, evolving in a potential W (whose typical shape is dis- 
played in Fig. [3]). The integration of Eq. ([9]) starts from 
the left, using the initial "position" A(X < 0) = 1, which 
corresponds to = W(l) = (supersonic uniform flow, 
see [16|). The behavior of A for X > L depends on the 



final value of A stationary solution exists only if 

A(X > L) remains bounded, i.e., < W(A±) (where 

A\ corresponds to the local maximum of W, see Fig. [3]). 
In this case, the transmission coefficient is [22j 



T 



1 



1 + (2k 2 £, 2 y 1 E t 



JVi 



(10) 



This formula is non-perturbative in the deep supersonic 
regime v 3> 1 ( i.e., it is not limited by the condition 
\5n\ <C no) and T, given by Eq. (|T0|) . may become very 
small (contrarily to Eq. ([7]) which is only valid for values 
of the transmission close to unity). 

In order to describe the density profile in the disordered 
region and to determine the value of E 1 , one can device 
a diffusion equation describing the evolution of E c \ under 
the effect of the kicks caused by each randomly placed 
impurity (cf. Fig. [3]). This allows -through (fTO)) - to 
compute the disorder average of the transmission and its 
probability distribution. For L L\ oc and k ni,^ 1 
one obtains 



(InT) = —L/L\ oc (k) , 



(11) 



where L\ oc {n) is given by Eq. 0, whereas the probability 
distribution reads 



PQnT) 



-Moc 



exp 



-^toc 

~4L 



L 

L\oc 



InT 



(12) 
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Excellent agreement of the analytical formulas (fTTj) and 
(fT2| with the numerical results is found, cf. Figs. &b) 
and[5{d). It is interesting to note that the exponential de- 
crease (fTT|) of the mean transmission and the log normal 
distribution (fT2")) are characteristic features of AL of lin- 
ear waves with a localization length L\ oc (see, e.g., [23j|). 
Hence, "typical" AL is predicted to occur in supersonic 
interacting BEC systems. Interaction induces a modifi- 
cation of the wave vector in L\ oc : the expression |7]) for 
the localization length coincides with the non-interacting 
one, but computed for an effective interaction-dependent 
wave vector k, instead of k = mV/h. The repulsive inter- 
action diminishes the available kinetic energy, and there- 
fore reduces the localization length with respect to the 
non-interacting case (since k < k). 

The present analysis also explains the onset of time de- 
pendence observed numerically in the supersonic region 
at fixed V and increasing L (see Fig. [TJ. Once the av- 
erage "classical energy" (E^) exceeds the value W(j4i) 
(see Fig. [3]), the number of stationary solutions decreases 
dramatically. This occurs for L > L*(k) with 



L*(k) =L 1oc (k) In (v 2 /8) /2 



(13) 



This expression fits very well the numerically determined 
frontier of the stationary supersonic domain (see Fig.[l}. 

At moderate supersonic velocities Eq. (fl~3|) shows that 
L* is roughly of the same order as L\ oc and one cannot 
reach the regime L\ oc <C L < L* where the log normal 
distribution (fT2")) is observed. This is the reason why 
Fig. [U(b) displays the distribution at relatively large ve- 
locity (v = 30, L* /L\ oc ~ 2.4). However, the exponential 
decay ([TTj) of the average transmission is observed for 
all stationary supersonic velocities, with L\ oc given by 
Eq. ©, as confirmed by Figs, da) and^b). 

The w-dependence of the two relevant length scales 
L\ oc and L* is shown in Fig. 1. Supersonic stationary so- 
lutions exist above a critical velocity v\ corresponding to 
L*(ki) = 0. One gets v\ = y/8 ~ 2.8, which is the onset 
of the ohmic regime. In contrast, AL exists only within 
the stationary regime for L\ oc < L < L* (cf. Fig. [1]). 
This is possible only if v > t>2, where L*{k2) — L\ oc {n2), 
i.e., Vi = ey/8 ~ 7.7. However, for v > V2, the difference 
of L\ oc (k) with its non- interacting counterpart L\ oc {k) is 
less than 2%. Hence interactions play no significant role 
in the AL regime. They will be relevant however in the 
low-velocity, small-disorder-length sector of the ohmic 
regime. At v — 3 for instance, they induce a modifica- 
tion of 11% of the localization length [cf. the comparison 
of interacting and non interacting results in Fig. [5Ja)] . 

In summary, we provide a complete picture of the quasi 
ID transport of a weakly interacting BEC through a dis- 
ordered potential. We find that interactions have several 
major and experimentally relevant consequences : (i) ex- 
istence of a subsonic SF regime, (ii) existence of non 
steady flows around v ~ 1, (iii) renormalization of the 
localization length in the stationary supersonic regime, 



and (iv) introduction of a maximum disorder length scale 
L* at which AL disappears and time dependence sets in. 
Explicit methods to characterize experimentally the dif- 
ferent regimes through, e.g., the heating rate, will be 
discussed elsewhere. 
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